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Using the time-dependent Ginzburg-Landau equations we study the propagation of planar fronts in 
superconductors, which would appear after a quench to zero applied magnetic field. Our numeri- 
cal solutions show that the fronts propagate at a unique speed which is controlled by the amount 
of magnetic flux trapped in the front. For small flux the speed can be determined from the lin- 
ear marginal stability hypothesis, while for large flux the speed may be calculated using matched 
asymptotic expansions. At a special point the order parameter and vector potential are dual, leading 
to an exact solution which is used as the starting point for a perturbative analysis. 
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There exists a wide class of problems in which a sys- 
tem, subjected to a sudden destabilizing change, re- 
sponds by forming fronts which propagate into the unsta- 
ble state. This phenomena occurs in models of popula- 
tion dynamics pulse propagation |4| , liquid crystals 
solidification P|, and tubular vesicles 0|; further ex- 
amples are discussed in Ref. ||^. The interesting issues 
common to all of these problems are whether the fronts 
approach a constant, unique speed, and if so, how the 
system selects this speed out of a manifold of possible 
speeds. The prototype is Fisher's equation 



dtu 



F{u), 



(1) 



where u > may be interpreted as a population density, 
and F{0) = F{1) = 0. As shown rigorously by Aron- 
son and Weinberger for a sufficiently localized initial 
condition the solutions of this equation evolve into fronts 
of the form u{x, t) = U{x — ct) which connect u = 1 at 
X — —oo to u ~ ai X — oo; the speed of the front 
satisfies 2-y/F'(0) < c < 2 sup^/ F{u)/u, so that for the 



special case F{i 



the selected speed is c = 2. 



There have been many attempts to generalize these re- 
sults to more complicated equations and physical sys- 
tems, including heuristic methods such as the marginal 
stability hypothesis (MSH) |^,@ and the structural sta- 
bility hypothesis [nilJl3], construction of exact solutions 
[^,|3|, variational methods [p^ , p^|jr^ ], and dynamical 
systems methods |p^ . 

In this paper we will study a closely related problem 
of front (or interface) propagation in superconductors. 
The problem which we have in mind is the following: be- 
gin with a bulk superconductor in an applied magnetic 
field equal to the critical field He, so that there is a sta- 
tionary, planar superconducting-normal interface which 
separates the normal and superconducting phases; then 
rapidly reduce the applied field to zero, so that the now 
unstable interface propagates toward the normal phase so 
as to expel any trapped magnetic flux, leaving the sam- 
ple in the Meissner phase. Assuming that the interface 



remains planar, what is its dynamics? In this paper we 
will show that under these conditions constant velocity 
fronts do propagate, at a unique speed which is controlled 
by the amount of magnetic flux which is trapped in the 
front jlj]. We calculate this speed in different param- 
eter regimes using the MSH, a perturbative calculation 
about an exact solution which we have discovered, and 
an asymptotic analysis valid for small speeds. Where ap- 
propriate, our analytic results are compared to extensive 
numerical solutions of the dynamic equations. 

To analyze the behavior of the superconducting- normal 
boundary, we use the one-dimensional time-dependent 
Ginzburg-Landau (TDGL) equations, which in dimen- 
sionless units llql are 



dtf = -^dlf -q'f + f~f 



adtq ^ dlq - f\, 



(2) 



(3) 



where / is the magnitude of the superconducting order 
parameter, q is the gauge- invariant vector potential (such 
that h ^ dxq IS the magnetic field), k is the Ginzburg- 
Landau parameter, and a is the dimensionless normal 
state conductivity (the ratio of the order parameter dif- 
fusion constant to the magnetic field diffusion constant). 
Notice that if the vector potential is zero, then Eq. (g) 
is exactly Fisher's equation, Eq. (0), with fronts which 
propagate at a speed c = 2/k in our units. Since lengths 
are measured in units of the penetration depth (typically 
of order 500 A in type-I superconductors) and time in 
units of the order parameter relaxation time (of order 
10~^-10~^'' s), the characteristic scale for speeds is of or- 
der 100 m/s [|l9|. For the propagating solutions which we 
are considering (i.e., after the field quench), f — 1, q ^ Q, 
and h = Q as X ^ —oo (the superconducting phase), and 
f = 0, q = Qoc, and /i = as X ^ oo (the normal phase). 
The physical meaning of Qoo is clear once we notice that 
the integrated magnetic field in the front (i.e., the to- 
tal magnetic fiux per unit length parallel to the front) is 
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J-oo '^(•^) ~ important 
control parameter for the front dynamics — the larger the 
trapped magnetic flux in the front the smaller the front 
speed. 

We have solved Eqs. (H) and (^) numerically for a wide 
range of k, a and Qoo, using the Crank-Nicholson method 
pO| . An initial configuration for the order parameter and 
magnetic field with the appropriate boundary conditions 
is established and then allowed to evolve in time. For 
our boundary conditions (in particular, /i = as x ^ 
oo), the front rapidly approaches a constant velocity. On 
an IBM RS 6000/370 approximately 5000 cpu minutes 
are needed to trace this evolution (allowing 300-500 time 
units to elapse usually brings us sufficiently close to a 
constant velocity solution). We analyze the profile of 
the order parameter and magnetic field for the constant 
velocity solutions (see Fig. |l] for a representative result) 
and determine the value of the front velocity as a function 
of Qoo, as shown in Fig. ^. The remainder of this paper 
is devoted to an analysis of the TDGL equations which 
will shed some light on these numerical results. 

To begin our analysis we will search for steady travel- 
ing wave solutions of the TDGL equations, of the form 
f{x,t) = F{X) and q{x,t) = Q{X) where X ^ x - ct 
with c > 0. Then the TDGL equations become 
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F" + cF' - Q'^F + F~F'^ 



0, 



(4) 
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FIG. 2. Numerical values (open squares) of the front speed 
as a function of Qaa for k = 1 and ct = 1. Also shown are 
the results of the linear marginal stability analysis (solid line) 
and matched asymptotic analysis (dashed line). 

and d), we find F{X) - exp(-A_A:), with 



(6) 



Since the magnitude of the order parameter F > 0, we 
require Im(A_ ) = in order to prevent the solution from 
oscillating, which implies 



Q" + acQ' - F'^Q = 0, 



(5) 



(the primes denote differentiation with respect to X). 
The order parameter connects i^ = latX = — ooto 
F = at A = oo, while the vector potential connects 
Q = at A" = — oo to Q = Qoo at A' = oo. In the spirit 
of the MSH, we first examine the linear stability of the 
leading edge of the front (X oo). Linearizing Eqs. (Q) 
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FIG. 1. Numerical solution for it — 1, a = 1, and Qtxi = 2, 
resulting in a front which moves to the right with speed 
c = 0.386. Shown are the order parameter (solid line), vector 
potential (dashed line), and field (dashed-dotted line). 



c > c' 



(Qoo < 1), 



(7) 



with c > for Qoo > 1. In addition to this lower bound 
on the speed, we conjecture that c < 2/k, since we expect 
the front to have the greatest speed in the absence of any 
trapped flux. The MSH |9l|lCt| would have us set c = c'; at 
this speed the leading edge of the order parameter front 
is marginally stable with respect to perturbations — the 
phase and group velocities of the perturbation are equal 
at this speed. At Qoo = 0, this yields c — 2/k, which is 
the rigorous result for Fisher's equation ||,^. In Fig.|| we 
compare the result of the MSH against the numerical re- 
sults. We see that for small Qoo there is a close agreement 
between the numerics and the MSH result. However, as 
Qoo approaches 1 we begin to see significant departures 
from the MSH; the MSH predicts that the front should 
be stationary for Qoo > 1. We therefore see that the 
linear MSH fails to accurately predict the selected front 
speed for sufficiently large trapped flux. 

Further understanding of the selection problem can be 
gleaned from an exact solution of the TDGL equations 
for Qoo — 1, K = l/\/2, and a — 1/2. To motivate this 
solution, we first notice that when Qoo = 1, it appears 
qualitatively that F + Q = 1. With this in mind, we 
look for solutions to (^) and (^) of the form F{X) = 
1 — Q{X); by substituting this Ansatz into (^) and (||), 
we find that this is only consistent when k — 1/ \/2 and 
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(T = 1/2. For this set of parameters the equation for the 
order parameter reduces to 



F" 



-F' + - = 0. 



(8) 



The exact solution of Eq. (||) can be constructed using 
the reduction of order method described by van Saarloos 
(see Ref. @, Sec. IV), with the result that c = V2, 
and the order parameter and vector potential are (up to 
translations) 



Fix) 



1 



Q{x) = 



1 



,-X/V2 ^ I ' 



(9) 



It should be emphasized that this duality between F and 
Q is distinct from the well known duality between the 
order parameter and the magnetic field which exists in 
the equilibrium Ginzburg-Landau (GL) equations at k = 
1/V2 m. To the best of our knowledge this is the only 
known exact solution of the TDGL equations. 

The exact solution serves as a useful check on our nu- 
merical results; our numerical solution of (||) and (||) for 
this set of parameters gives c — 1.4142, bolstering con- 
fidence in the numerical work. The exact solution also 
illustrates the limitations of the linear MSH, which pre- 
dicts a speed of c = for these parameters. Finally, 
the exact solution can also serve as the starting point for 
perturbative solutions of the TDGL equations, along the 
lines of the renormalization group method discussed in 
Ref. [O. We find I2I 



= V2 



1 



1 2 



(10) 



with = 2 - I/k^, Ecr = c - 1/2, and eg = Qoo - 1- 
For example, when k — 0.6, a — 0.5, and Qoo = 1, 
Eq. ( p^ ) yields c — 1.5975, while the numerical result is 
c= 1.5826. 

Although the exact solution described above is a use- 
ful touchstone for our numerical work, we desire a more 
general approach to the problem of calculating the front 
speed. To do this we develop a perturbative solution of 
Eqs. (^) and (||) for small c (large trapped flux); a simi- 
lar treatment for curved interfaces in two dimensions has 
been given in Refs. ||l8| and Expand the solutions 

in powers of c (the inner expansion) , 



F{X;c) = Fo{X) + cF,{X) + 
Q{X;c)=Qo{X) + cQi{X) + 



(11) 



and substitute these expansions into Eqs. ^ and (||). 
The 0(1) equations for (Fq, Qq) are the GL equations. 
The important feature of the solutions for our purposes 
is that the vector potential in the normal phase {X —^ 00) 
has the asymptotic behavior (up to a translation) 



e.s.t.. 



(12) 



where e.s.t.= "exponentially small terms." This result 
shows that the magnetic field in the normal phase ap- 
proaches l/-\/2 (which is He in conventional units) — a 
planar interface can only be in equilibrium when the field 
in the normal phase is He- Proceeding to 0(c), we have 

+ [l-3(Fo)2]Fi = -F^, (13) 



Q'l - (Fo)2Qi - 2F0Q0F1 = -aQ'o. (14) 

The asymptotic behavior of Qi can be obtained as fol- 
lows. Multiply Eq. (|l|) by F(^, Eq. (|l|) by Q'^, add the 
two equations together, and integrate the result from —00 
to X. Then integrate F{' and Q'{ by parts twice, and use 
the fact that (Fq^Qq) are solutions to the homogeneous 
versions of Eqs. ([l3| ) and (|^ (the zero mode). The final 
result for the asymptotic behavior as X — + cx3 is 



X + Ci+ e.s.t. 



2^/2 x/2 
where Ci is an integration constant and /? is given by 



(15) 



/3 = 2 



iF[,)^ + aiQ',r-^Q', 



dx. 



(16) 



The kinetic coefficient /3(k, a) must be determined nu- 
merically from the solutions of the GL equations . 

By comparing Qo and cQi, we see that these two terms 
become comparable when cX = 0(1), indicating a break- 
down of the perturbative expansion. This suggests intro- 
ducing the owter variable ^ = cX, with /(^) = F{X) and 
liO — Q{X)- In terms of these outer variables Eqs. (^) 
and M) become 



c'q" 



^c\' - fq - 0. 



(17) 



(18) 



Expanding the solutions in powers of c, 

/(e;c) = /o(0 + c/i(0 + --- 
g(C;c) = h gi(0 + • ■ • 



(19) 



we find that /o = 1, /i = in the outer superconducting 
region (^ ^ — 00), /o = /i = in the outer normal region 
(^ — > 00), and go = 91 = in the outer superconducting 
region, with 

90(0 = Ao + Soe-*«, gi(C) = Ai+Bie-*« (20) 

in the outer normal regions; An and i?„ are integra- 
tion constants. The inner and outer solutions must now 
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be matched together in an appropriate overlap region 
M-ESi, with the result 



= -Bo 



1 



/3 



(21) 



Using this expansion, we can determine the asymptotic 
value of the vector potential in the normal phase as an 
expansion in powers of c: 



1 



p_ 



0(c). 



If we think of Qoo as the control variable, then we have 
for the selected velocity 



1 



(23) 



As was pointed out in the kinetic coeffi- 

cient may actually be negative for some values of the 
parameters. One might worry that when — —P^ 

Eq. ( p3|) predicts an infinite velocity, contradicting our 
conjecture that c < 2//t. However, if this occurs we sim- 
ply have a breakdown of the perturbation expansion, in- 
dicating the need to keep higher order terms. In Fig. |^ we 
compare the asymptotic result, Eq. (p3|), and the numer- 
ical results, for k = \ and ct = 1; the kinetic coefficient 
for these parameters is /3 = —0.216 24 1. The agreement 
is excellent in the appropriate region of large Qoo ■ The 
agreement is equally impressive for other values of k and 
a. In particular, when k — l/\/2 and a = 1/2, /3 = [ p^ , 
so that the asymptotic analysis predicts c = V^/Qoo', 
if we set Qoo = 1 + cq, and expand for small eg, then 
c = -\/2[l — eQ + 0(e^)], which agrees to lowest order with 
the expansion about the exact result given in Eq. (|lO|). 

In summary, we have studied the propagation of fronts 
separating the superconducting and normal phases, 
which are produced after a quench to zero applied mag- 
netic field. In addition to its possible practical impor- 
tance in understanding flux expulsion in superconduc- 
tors, this problem provides an interesting variation on 
the theme of front propagation in unstable systems. By 
varying the amount of flux trapped in the front we can 
go continuously from a regime at small flux in which the 
speed is close to that predicted by the linear marginal 
stability hypothesis, to a regime at high flux which can 
be treated using the method of matched asymptotic ex- 
pansions. While we have no general results at intermedi- 
ate values of the flux, for a particular set of parameters 
we have discovered an exact solution, which serves as 
the starting point for a perturbative calculation in this 
regime. We are currently expanding our study to in- 
clude front propagation in two dimensions and a detailed 
study of the diffusive fronts which appear after quenches 
to non-zero applied flelds p2[ . 
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